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The Sandage-Loeb (SL) test is a unique method to explore dark energy at the "redshift desert" 
(2 < z < 5), an era not covered by any other dark energy probes, by directly measuring the temporal 
variation of the redshift of quasar (QSO) Lyman-a absorption lines. In this paper, we study the 
prospects for constraining the new agegraphic dark energy (NADE) model and the Ricci dark energy 
(RDE) model with the SL test. We show that, assuming only a ten-year survey, the SL test can 
constrain these two models with high significance. 

I. INTRODUCTION 

The astronomical observations of type la supernovae (SNIa) indicate that our universe is undergoing an 
accelerating expansion [1]. This cosmic acceleration has also been confirmed by other observations, such 
as the large scale structure (LSS) [2] and the cosmic microwave background (CMB) [3]. Nowadays it is 
the most accepted idea that a mysterious dominant component, dark energy, with large enough negative 
pressure, is responsible for this cosmic acceleration. Among all theoretical models, the preferred one is the 
so-called ACDM model, which consists of a mixture of Einstein's cosmological constant A and the cold 
dark matter (CDM). The ACDM model provides an excellent explanation for the acceleration of the universe 
and the existing observational data. However, the cosmological constant has to face severe theoretical 
problems such as the puzzle why the dark energy density today is so small compared to typical particle 
scales. Therefore, except the ACDM model, many dynamical dark energy models have been proposed, in 
which the equation of state (EoS) of dark energy is no longer a constant but slightly evolves with time. For 
reviews of dark energy, see, e.g., Ref. [4]. 

In the face of so many candidate models, it is extremely important to identify which one is the correct 
model by using the observational data. The measurement of the expansion rate of the universe at diff"erent 
redshifts is crucial to discriminate these competing candidate models. Up to now, a number of cosmological 
tools have been used to successfully probe the expansion and the geometry of the universe. These, typically, 
include the luminosity distance of SNIa, the position of acoustic peaks in the CMB angular power spectrum, 
and the scale of baryon acoustic oscillations (BAO) in the power spectrum of matter extracted from galaxy 
catalogues. Recently, the application of time evolution of cosmological redshift as a test of dark energy 
models has become attractive, since this method opens a new window of exploring the "redshift desert" 
(2 < z < 5). In addition to being a direct probe of the dynamics of the expansion, this method has the 
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advantage of not relying on a determination of the absolute luminosity of the observed sources, but only 
on the identification of stable spectral lines, so this method can reduce the uncertainties from systematic or 
evolutionary effects. 

Sandage [5] was the first to propose the possible application of this kind of observation as a cosmological 
tool. However, owing to the tininess of the expected variation, this observation was deemed impossible at 
that time. In 1998, Loeb [6] revisited this suggestion and argued that the redshift variation of quasar (QSO) 
Lyman-a absorption lines could be detected in the not too distant future, given the advancement in technol- 
ogy occurred over the last forty years. In fact, the cosmological redshift variation at Icr would be detected 
in a few decades, if a sample of a few hundred QSOs could be observed with high resolution spectroscopy 
with a ten meter telescope. This method is usually refeiTcd to as "Sandage-Loeb" (SL) test. The possibil- 
ity of detecting the temporal variation of redshift with the Cosmic Dynamics Experiment (CODEX) was 
first analyzed by Corasaniti, Huterer and Melchiorri [7]. Their work [7] has provided the first quantitative 
analysis of the SL test, from which all other analyses have followed. 

In Ref. [7], Corasaniti, Huterer and Melchiorri employed the SL test to constrain dark energy models 
such as ACDM model, Chaplygin gas model, and interacting dark energy model. Later, Balbi and Quer- 
cellini [8] extended this analysis to more dark energy models including constant EoS model, variable EoS 
model, interacting dark energy model, DGP model, Cardassian model, generalized Chaplygin gas model, 
affine EoS model, etc. More recently, Zhang, Zhong, Zhu and He [9] further used the SL test to explore 
the holographic dark energy model. However, it should be pointed out that there are three holographic dark 
energy models: the original holographic dark energy model [10], the new agegraphic dark energy model 
[11], and the Ricci dai^k energy model [12]. Actually, in Ref. [9], only the original holographic dark energy 
model [10], i.e., the model in which the IR cutoff is given by the future event horizon, was investigated. 
Thus, along this line, as a next step, one should further explore the new agegraphic and the Ricci dark energy 
models with the SL test. In this paper, this will be done. This will provide a complementary to the work of 
Ref. [9] and keep the investigation of holographic dark energy models more complete. 

In the subsequent section, we will briefly review the new agegraphic dark energy model and the Ricci 
dark energy model. In Sec. Ill, we will explore these two models with the SL test. In the last section, we 
will give some concluding remarks. 

II. NEW AGEGRAPHIC AND RICCI DARK ENERGY MODELS 

In this section, we will briefly review the new agegraphic dark energy model and the Ricci dark energy 
model. In fact, these two models both belong to the holographic scenario of dai^k energy. 
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It is well known that the holographic principle is an important result of the recent research for exploring 
the quantum gravity [13]. This principle is enlightened by investigations of the quantum property of black 
holes. In a quantum gravity system, the conventional local quantum field theory will break down because 
it contains too many degrees of freedom that would lead to the formation of a black hole breaking the 
effectiveness of the quantum field theory. To reconcile this breakdown with the success of local quantum 
field theory in describing observed particle phenomenology, some authors proposed a relationship between 
the ultraviolet (UV) and the infrared (IR) cutoffs due to the limit set by the formation of a black hole [14]. 
The UV-IR relation in turn provides an upper bound on the zero-point energy density. In other words, if the 
quantum zero-point energy density pvac is relevant to a UV cutoff, the total energy of the whole system with 
size L should not exceed the mass of a black hole of the same size, and thus we have L^Pvac ^ LMp^ This 
means that the maximum entropy is of the order of S^^. When we take the whole universe into account, the 
vacuum energy related to this holographic principle is viewed as dark energy, usually dubbed holographic 
dark energy (its density is denoted as pde hereafter). 

The largest IR cutoff L is chosen by saturating the inequality so that we get the holographic dark energy 
density [10] 

Pde = 3c2m2jL-2 , (1) 

where c is a numerical constant, and Mpi = 1/ VSttG is the reduced Planck mass. If we take L as the size of 
the current universe, for instance the Hubble radius H~^, then the dark energy density will be close to the 
observational result. However, if one takes the Hubble scale as the IR cutoff, the holographic dark energy 
seems not to be capable of leading to an accelerating universe [15]. The first viable version of holographic 
dark energy model was proposed by Li [10]. In this model, the IR length scale is taken as the event horizon 
of the universe. The holographic dai^k energy model based on the event horizon as the IR cutoff has been 
widely studied [16] and found to be consistent with the observational data [17, 18]. 

There are also other two versions of the holographic dark energy, i.e., the new agegraphic dark energy 
model [11, 19] and the Ricci dark energy model [12, 20, 21]. For the new agegraphic dark energy model, 
the IR scale cutoff is chosen to be the conformal age of the universe; for the Ricci dark energy model, the IR 
cutoff is taken as the average radius of the Ricci scalar curvature. We shall briefly review these two models 
in the following subsections. 
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A. New agegraphic dark energy model 

For a spatially flat (the assumption of flatness is motivated by the inflation scenario) Friedmann- 
Robertson- Walker (FRW) universe with matter component pm and dark energy component pde, the Fried- 
mann equation reads 

3MI^H^ =p,^+Pi,, (2) 



or equivalently, 



where H = a/ais the Hubble parameter, 0.^0 is the present fractional matter density, and Ode = tt ~ -^m^''m2 
is the fractional dark energy density. 

In the old version of the agegraphic dark energy model [22], the IR cutoff is chosen as the age of the 
universe T (here it should be pointed out that the light speed has already been taken as 1, so time and length 
have the same dimension). However, there are some inner inconsistencies in this model; for details see 
Ref. [11]. So, in this paper, we only discuss the new version of the agegraphic dark energy model. In the 
new agegraphic dark energy model, the IR cutofi" is chosen to be the conformal age of the universe, 

Jo a Jo a^H ' 
so the density of the new agegraphic dark energy is 

Pde - 3«2m2 77-2. (5) 

To distinguish from the original holographic dark energy model, a new constant parameter n is used to 
replace the former parameter c. Taking derivative for Eq. (5) with respect to x = \na and making use of Eq. 
(4), we get 



Pde ^ -2pde . (6) 

na 



This means that the EoS of the new agegraphic dark energy is 



2 V^dT 

Wde = -1 + :; ■ (7) 

5n a 



Taking derivative for Q.dc = n^/iH^rj^), and considering Eq. (4), we obtain 



fi^, = 2ade f--^ ' (8) 
na 



5 



where 



3 3 3 

e=-{l+ WdeHde) - - - -Ode + — - (9) 

2 2 2 na 



Hence, we get the equation of motion for Qde> 

2 Vf^de 



i^de = ^^de(l - i^de) 3 —\, (10) 



and this equation can be rewritten as 



= -OdeCl -Hde) 3(1 +Z) 



|3(1+z)-1-^Vq^J . (11) 



dz 

As in Ref. [11], we choose the initial condition, Qde(Zini) - n^i^ +Zim)~^/4, at Zini = 2000, then Eq. (11) can 
be numerically solved. Substituting the results of Eq. (11) into Eq. (3), the function E(z) can be obtained. 
Notice that once n is given, QmO - 1 - ^deiz - 0) can be natural obtained by solving Eq.(ll), so the new 
agegraphic dark energy model is a single-parameter model. 



B. Ricci dark energy model 

For a spatially flat FRW universe, the Ricci scalar- is 

'R = -6(h + 2H^). (12) 

As suggested by Gao et al. [12], the energy density of Ricci dark energy is 

p,^ = i^(// + 2//2) = (13) 

where a is a positive numerical constant to be determined by observations. Comparing to Eq. (1), it is seen 
that if we identify the IR cutoff LT^ with -"R/G, we have a = c^. As pointed out by Cai et al. [23], the 
Ricci dark energy can be viewed as originated from taking the causal connection scale as the IR cutoff" in the 
holographic setting. Now, the Friedmann equation, in a universe containing Ricci dark energy and matter, 
can be written as 

„^- = ?|ip„„,-3.,,(lf?|!,2„.). (14) 

and this equation can be further rewritten as 

E' = 0^oe-'' + ai^^^ + 2E'y (15) 

where E = H/Hq. Solving this equation, and using the initial condition £0 = E{to) = 1, we have 

E(z) = (^(1 + z? + (1 - ^)(1 + zr~'A'" . (16) 
\2 - a 2 - a j 

There are two model parameters, £1^^ and a, in the Ricci dark energy model. 
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III. THE SANDAGE-LOEB TEST 

In this section, we will first review the Sandage-Loeb test, and then explore the new agegraphic dark 
energy model and the Ricci dark energy model with the SL test. 

First, let us consider an isotropic source emitting at rest. The well-known redshift relation of the radiation 
emitted by the source at ts and observed at to is 

aits) 

Furthermore, consider lights emitted after a period A?., at tg + Atg and detected later at t„ + Ato- Obviously, 
the observed redshift of the source at to + Atg is 

a(to + Ato) 

Zs{to + Ato)^ )' .; -1. (18) 
a{ts + Ats) 

Therefore, the variation of the source redshift between times to and to + Atg would be measured as follows: 

a{to + Ato) a{to) .. 

Azs = — — -— ■ (19) 

a{ts + Atg) a{ts) 

We can expand the ratio a(to + Ato)la(ts + A?s) to linear- order, under the approximation Atjt «; 1. Further- 
more, using the the relation Ato - \o.(to)la(ts)\Ats, we obtain 

a(to) - a{ts) 



Ato. (20) 



a{ts) 

It shows that the redshift variation Azs is directly related to a change in the expansion rate during the 
evolution of the universe, and it is thus a direct probe of the dynamics of the cosmic expansion. This 
redshift variation can be related to a spectroscopic velocity shift, Av = Azv/(1 + Zs)- Using the Hubble 
parameter H{z) = a{z)/a{z), we obtain 

E(Zs) 



Av = HoAto 



1 



(21) 



1 

where Hq is the Hubble constant and E{z) = H{z)/Hq. The function E{z) contains all the details of the cos- 
mological model under investigation. It is clear that the expansion history E{z) is related to the spectroscopic 
velocity shift via Eq. (21). 

Though the amplitude of the velocity shift is very small, the absorption lines in the quasar Lyman-o- pro- 
vide us with a powerful tool to detect such a small signal. Monte Cai^lo simulations of Lyman-or absorption 
lines have been performed to estimate the uncertainty on Av as measured by the CODEX spectrograph [24]. 
The statistical eiTor can be estimated as 



/2350\ 1 30 / 5 cm 
[s/nJ \N^[l+ZQsol 



^Av-1.4| — M— 1^ I — , (22) 
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where S IN denotes the spectral signal-to-noise defined per 0.0125 A pixel, Nqso is the number of Lyman-a 
quasars, and zqs o is the quasar's redshift. In order to detect the cosmic signal, a large S /N is necessary, 
but this implies that a positive detection is not feasible with current telescopes. Fortunately, the CODEX 
under design will be installed on the ESO Extremely Large Telescope. The necessary signal-to-noise can 
be reached by such an about 50 meter giant with just few hours integration. The velocity shift measure- 
ments open a cosmological window with particular- focus on dai^k energy models. From the velocity shift 
measurements, one can forecast constraints on parameters of cosmological models. In this paper, follow- 
ing Refs. [7-9], we consider experimental configuration and uncertainties similar to those expected from 
CODEX. We assume that the survey would observe a total of 240 QSOs uniformly distributed in six equally 
spaced redshift bins in the range 2 < z < 5, with a signal-to-noise S /N = 3000, and the expected uncertainty 
as given by Eq. (22). Also, in this paper, we consider a ten-year survey, namely, Ato = 10 years. In what 
follows, we shall explore the new agegraphic dark energy model and the Ricci dark energy model with the 
SL test. 




First, we discuss the velocity shift behavior in the new agegraphic dark energy (NADE) model. Numer- 
ically solving the differential equation (11) and inserting Eq. (3) into Eq. (21), one can reconstruct the Av(z) 
curves for the NADE model. Note that the NADE model is a single-parameter model, and the sole model 
parameter is n. The current cosmological constraint on the parameter nh n = 2.807^^°^^ (Icr) ^gj^Q (2o") 
[17]. In this work, we take the values of n as the central value as well as the one-sigma and two-sigma 
limits of the above observational constraint. In Fig. 1, we reconstruct the velocity shift behavior of the 
NADE model. In the left panel, we use the ACDM model as a fiducial model to perform an SL test. For 
the ACDM cosmology, Qmo is chosen to be 0.274 as given by the WMAP five-year observations [25]. Note 
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FIG. 2: The SL test for the RDE model. The ACDM model with = 0.274 is used to be the fiducial model. In the 
left panel, we fix Q.^q - 0.324 and vary a; in the right panel, we fix a = 0.37 1 and vary Q^o- 
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FIG. 3: The SL test for the RDE model. The RDE model with il^o = 0.324 and a = 0.371 is used to be the fiducial 
model. In the left panel, we fix Qmo - 0.324 and vary a; in the right panel, we fix a = 0.371 and vary Q.mO- 



also that the Hubble constant Hq in Eq. (21) is taken as 72 km s~^ Mpc~\ in the whole discussion. From 
the left panel of Fig. 1 we see that the NADE model can be distinguished from the ACDM model via the 
SL test. However, it is interesting to notice that the curve with n = 0.2894 still lies in the Icr error range 
of the ACDM fiducial model within the redshift desert 2 < z < 5. So, it is fair to say that actually the SL 
test could not completely distinguish the NADE model from the ACDM model, though it can do it rather 
effectively. Of course, if the SL test is combined with the low-redshift observations such as the SNIa, weak 
lensing, and BAO, the NADE model can be completely distinguished from the ACDM model. On the other 
hand, we also perform an SL test by using the NADE model with n = 2.807 (the central value of n given by 
the current observational constraint) as a fiducial model. We show this case in the right panel of Fig. 1. In 
this case, we see that the prospective SL test is very powerful to be used to constrain the NADE model, and 
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it is cleai^ly better than the current low-redshift observations. 

Next, we switch to the Ricci dark energy (RDE) model. For the RDE model, the Hubble expansion 
history is described by Eq. (16), so the Av(z) behavior can be reconstructed by substituting Eq. (16) into 
Eq. (21). The RDE model has two model parameters, a and Q.mO, and the current observational constraint 
results are [17]: a = 0.37l!°°23 (j^) +0.037 ^2cr) and Q^q = 0.324+0024 (^^^ +o.040 ^2cr). Like the above 

discussion, we will still first employ the standard dark energy cosmology, the ACDM model, as a fiducial 
model to perform an SL test for the RDE model. For the fiducial ACDM model, we still take QmO = 0.274, 
the value given by WMAP [25]. In this test, one can forecast how well a deviation of the RDE model from 
the ACDM model can be detected. We show the results of such an SL test in Fig. 2. In the left panel, the 
parameter 0^0 is fixed to be 0.324 for the RDE model, and the parameter a of the RDE model is adjustable; 
in the right panel, we fix o- = 0.371 and vary Q.^q for the RDE model. From this figure, it is clear to 
see that the SL test in the "redshift desert" is very successful in distinguishing the RDE model from the 
standard ACDM cosmology. The two models can be distinguished via the SL test at about 6 - 8 o" level, 
assuming only a 10-year survey. This result is reasonable and understandable, since the RDE model has a 
tracker behavior [20] so that its matter-dominated era is different from that of the ACDM model, and the 
SL test is just the best way to probe the behavior of the matter-dominated phase of dark energy models. 
Furthermore, it is worth noticing that such measurements are mostly sensitive to the matter density f2niO> 
while the dependence on a is weaker. To forecast whether the SL test is able to break the degeneracy of 
the parameters of the RDE model, we further use the RDE model with QmO - 0.324 and a = 0.371 (the 
central values of the current observational constraint) as the fiducial model to perform an SL test. The 
results of this test are plotted in Fig. 3. In the left panel, we fix 0,^,0 = 0.324 and vary a in the 2 cr range 
of the current observational result; in the right panel, we fix a = 0.371 and vary QmO in the 2 cr range of 
the current observational result. From Fig. 3, we can see that the SL test with a 10-year survey could not 
provide precision determination of the parameter a, but could precisely determine the value of Q^q for the 
RDE model. This result is in agreement with the previous work on the SL test [7]; in Ref. [7] it is shown 
that the Sandage-Loeb constraints on w are not competitive with those inferred from other low-redshift 
observations. The equation of state w of RDE is mainly determined by the parameter a, so the SL test 
could not constrain the value of a precisely. However, as pointed out in Ref. [7], one should note that the 
constraints obtained by SL test decrease linearly with time, so for measurements made over a century, and 
with the expected larger number of QSOs, the SL limits on w can easily be at the few percent level. 

Finally, let us show how well the proposed SL test could constrain the two dark energy models by using 
a Fisher matrix method. To find the expected precision of the SL test with CODEX, one must assume a 
fiducial model, and then simulate the experiment assuming it as a reference model. We employ the ACDM 



10 




FIG. 4: Predicted likelihood distribution of the parameter n of the NADE model from the SL test. 
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FIG. 5: Predicted probability contours at 68.3% and 95.4% confidence levels in the (Q.,„o, Q") plane for the RDE model 
from the SL test. 



model as the fiducial model and produce the mock data of the velocity-drift in the redshift desert with the 
error bars given by Eq. (22). Next, we will perform a Fisher matrix analysis on the model parameter space 
with the above assumption. The of the analysis is given by 



240 



.data/ 



(23) 



where Av''^'^(z,) denotes the mock data produced by the fiducial model, Av™°'^^\zi) is the theoretical predic- 
tion of the dark energy model under investigation, and o"^^,(Z() is the eiTor bar estimated from Eq. (22). The 
NADE model is a single-parameter model, so once the parameter n is given, the parameter Qmo as well as 
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other properties and dynamical evolution will be determined accordingly. In Fig. 4 we plot the likelihood 
distribution of the parameter n of the NADE model as expected from the SL test. We see that the parameter 
n is able to be accurately determined by the SL test, n = ^■'^54^^-^'^ (Ic)- So, the SL test could obtain 
cr„ 0.005, better than the current observational result by at least one order of magnitude. The RDE model 
has two model parameters, Q^o and a. In Fig. 5 we plot the Icr and Icr contours in the QmO - en plane for the 
RDE model. The Icr results are: Q^o - 0.2055+|J °°^^ and a = 0.5000+° {]262 ^^le SL test would obtain 
Cfimo - 0.005, better than the cun^ent constraint by one order of magnitude, and ctq, - 0.03, at the same 
accuracy level compared to the current observational result. The above analysis reinforces the conclusion 
that the two dark energy models indeed can be constrained by the SL measurements with high significance. 



rv. CONCLUDING REMARKS 



The Sandage-Loeb test is a promising method for constraining dark energy by using the direct measure- 
ments of the temporal shift of the quasar Lyman-or absorption lines at high redshift (2 < z < 5). While 
the signal of this effect is extremely small, the near-future large telescopes with ultrahigh resolution spec- 
trographs (such as the CODEX under design) will definitely be capable of measuring such a small signal 
over a period as short as ten years. Notwithstanding, one still may ask whether it is worthwhile to probe 
z > 2, since in the standard dark energy cosmology such as the ACDM model or the slowly rolling scalar 
field model the dark energy is subdominant at redshift z > 1 and almost completely negligible at z > 2. The 
answer to this question is affirmative. As explained in Ref. [7], it is quite rational to look for the signatures 
of dark energy at all available epochs, because at present we do not know much about the physical nature 
and cosmological origin of dark energy, and there are too many possibilities for dark energy. Especially, 
there are lots of dark energy models in which dark energy density is non-negligible at high redshift. So, 
it is meaningful to study the future observations of velocity shift and their impact on dark energy models. 
The SL test on Chaplygin gas model and interacting dark energy model has been studied in Ref. [7]. In 
Ref. [8], more classes of dark energy models have been investigated with the SL test. In particular, the 
original holographic dai^k energy model has been explored with the SL test in Ref. [9]. So, the exploration 
of the other two holographic dark energy models, the new agegraphic dark energy model and the Ricci dark 
energy model, is naturally the next step, and this will provide a complementary to the work of Ref. [9] and 
keep the investigation of holographic scenarios of dark energy more complete. 

In this paper, we have analyzed the prospects for constraining the new agegraphic dark energy model 
and the Ricci dark energy model at the redshift desert 2 < z < 5 from the Sandage-Loeb test. The NADE 
model is a sole-parameter model; the parameter n together with an initial condition can determine the whole 
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cosmological evolution history of the NADE model. Actually, the evolution of this model is similar to that 
of a slowly rolling scalar field. Thus, though the SL test can be used to distinguish the NADE model from 
the ACDM model, this discrimination is not absolute. It is obvious that if the SL test is combined with the 
low-redshift observations such as the SNIa, weak lensing, and BAO, the NADE model can be completely 
distinguished from the ACDM model. Furthermore, we forecast that the prospective SL test is very powerful 
to constrain the NADE model, and it is better than the cuiTcnt low-redshift observations. The RDE model 
is a two-parameter model, and it has a tracking solution so that it is very different from the ACDM model 
in the matter-dominated epoch. So, the SL test is very successful in distinguishing the RDE model from the 
ACDM model. The two models can be distinguished via the SL test at about 6 - 8 o" level, assuming only 
a 10-year survey. As the SL test mostly probes the matter density at high redshift, the constraint on 0^0 is 
very strong, while the constraint on a is much weaker. 

Acknowledgments 

We are grateful to the referee for the insightful comments and suggestions, which have led to significant 
improvement to this paper. We also thank Dr. Hongbao Zhang for useful discussions. This work was 
supported by the Natural Science Foundation of China under Grant Nos. 10705041 and 10975032. 



[1] A. G. Riess et al. [Supernova Search Team Collaboration], Astron. J. 116, 1009 (1998) [arXiv:astro- 
ph/9805201]; S. Perlmutter et al. [Supernova Cosmology Project Collaboration], Astrophys. J. 517, 565 (1999) 
[arXiv:astro-ph/9812133]. 

[2] M. Tegmark et al. [SDSS Collaboration], Phys. Rev. D 69, 103501 (2004) [arXiv:astro-ph/0310723]; K. Abaza- 
jian et al. [SDSS Collaboration], Astron. J. 128, 502 (2004) [arXiv:astro-ph/0403325]. 

[3] D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 148, 175 (2003) [arXiv:astro-ph/0302209]; 
C. L. Bennett efflZ. [WMAP Collaboration], Astrophys. J. Suppl. 148, 1 (2003) [arXiv:astro-ph/0302207]. 

[4] V. Sahni and A. A. Starobinsky, Int. I. Mod. Phys. D 9 (2000) 373 [arXiv:astro-ph/9904398]; S. M. Carroll, 
Living Rev. Rel. 4 (2001) 1 [arXiv:astro-ph/0004075]; P J. E. Peebles and B. Ratra, Rev. Mod. Phys. 75 
(2003) 559 [arXiv:astro-ph/0207347]; T Padmanabhan, Phys. Rept. 380 (2003) 235 [arXiv:hep-th/0212290]; 
E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D 15 (2006) 1753 [arXiv:hep-th/0603057]; 
E. V. Linder, Am. J. Phys. 76, 197 (2008) [arXiv:0705.4102 [astro-ph]]. 

[5] A. Sandage, Astrophys. J. 139, 319 (1962). 

[6] A. Loeb, Astrophys. J. 499, LI 1 1 (1998) [arXiv:astro-ph/9802122]. 

[7] P S. Corasaniti, D. Huterer and A. Melchiorri, Phys. Rev. D 75, 062001 (2007) [arXiv:astro-ph/0701433]. 
[8] A. Balbi and C. Quercellini, Mon. Not. Roy Astro. Soc. 382, 1623 (2007) [arXiv:0704.2350 [astro-ph]]. 



13 



[9] H. Zhang, W. Zhong, Z. H. Zhu and S. He, Phys. Rev. D 76, 123508 (2007) [arXiv:0705.4409 [astro-ph]]. 
[10] M. Li, Phys. Lett. B 603, 1 (2004) [arXiv:hep-th/0403127]. 

[11] H. Wei and R. G. Cai, Phys. Lett. B 660, 113 (2008) [arXiv:0708.0884 [astro-ph]]. 

[12] C. Gao, F. Q. Wu, X. Chen and Y. G. Shen, Phys. Rev. D 79, 04351 1 (2009) [arXiv:0712.1394 [astro-ph]]. 
[13] G. 't Hooft, arXiv:gr-qc/93 10026; L. Susskind, J. Math. Phys. 36, 6377 (1995) [arXiv:hep-th/9409089]. 
[14] A. G. Cohen, D. B. Kaplan and A. E. Nelson, Phys. Rev. Lett. 82, 4971 (1999) [arXiv:hep-th/9803132]. 
[15] S. D. H. Hsu, Phys. Lett. B 594, 13 (2004) [arXiv:hep-th/0403052]. 

[16] Q. G. Huang and M. Li, JCAP 0408, 013 (2004) [arXiv:astro-ph/0404229]; Q. G. Huang and M. Li, JCAP 0503, 
001 (2005) [arXiv:hep-th/0410095];X. Zhang, Int. J. Mod. Phys. D 14, 1597 (2005) [arXiv:astro-ph/0504586]; 
X. Zhang, Phys. Lett. B 648, 1 (2007) [arXiv:astro-ph/0604484]; X. Zhang, Phys. Rev. D 74, 103505 (2006) 
[arXiv:astro-ph/0609699]; J. Zhang, X. Zhang and H. Liu, Phys. Lett. B 651, 84 (2007) [arXiv:0706.1 185 
[astro-ph]]; Y. Z. Ma and X. Zhang, Phys. Lett. B 661, 239 (2008) [arXiv;0709.1517 [astro-ph]]; X. Wu and 
Z. H. Zhu, Phys. Lett. B 660, 293 (2008) [arXiv:07 12.3603 [astro-ph]]; X. Zhang, Phys. Lett. B 683, 81 (2010) 
[arXiv:0909.4940 [gr-qc]]; M. R. Setare, J. Zhang andX. Zhang, JCAP 0703, 007 (2007) [arXiv;gr-qc/061 1084]; 
J. Zhang, X. Zhang and H. Liu, Eur. Phys. J. C 52, 693 (2007) [arXiv:0708.3121 [hep-th]]; J. Zhang, X. Zhang 
and H. Liu, Phys. Lett. B 659, 26 (2008) [arXiv:0705.4145 [astro-ph]]; Y. G. Gong, Phys. Rev. D 70, 064029 
(2004) [arXiv:hep-th/0404030]; B. Wang, E. Abdalla and R. K. Su, Phys. Lett. B 611, 21 (2005) [arXiv:hep- 
th/0404057]; X. Wu, R. G. Cai and Z. H. Zhu, Phys. Rev. D 77, 043502 (2008) [arXiv:0712.3604 [astro-ph]]; 
M. Li, C. Lin and Y. Wang, JCAP 0805, 023 (2008) [arXiv:0801.1407 [astro-ph]]; M. Li, X. D. Li, C. Lin 
and Y. Wang, Commun. Theor Phys. 51, 181 (2009) [arXiv:081 1.3332 [hep-th]]; B. Wang, Y. G. Gong and 
E. Abdalla, Phys. Lett. B 624, 141 (2005) [arXiv;hep-th/0506069]; B. Wang, C. Y. Lin and E. Abdalla, Phys. 
Lett. B 637, 357 (2006) [arXiv:hep-th/0509107]; B. Wang, C. Y. Lin, D. Pavon and E. Abdalla, Phys. Lett. 
B 662, 1 (2008) [arXiv:07 11.2214 [hep-th]]; S. Nojiri and S. D. Odintsov, Gen. Rel. Grav. 38, 1285 (2006) 
[arXiv:hep-th/0506212]. 

[17] M. Li, X. D. Li, S. Wang and X. Zhang, JCAP 0906, 036 (2009) [arXiv:0904.0928 [astro-ph.CO]]. 

[18] X. Zhang and E Q. Wu, Phys. Rev. D 72, 043524 (2005) [arXiv:astro-ph/0506310]; X. Zhang and R Q. Wu, 
Phys. Rev. D 76, 023502 (2007) [arXiv:astro-ph/0701405]; Q. G. Huang and Y. G. Gong, JCAP 0408, 006 
(2004) [arXiv:astro-ph/0403590];Z. Chang, E Q. Wu and X. Zhang, Phys. Lett. B 633, 14 (2006) [arXiv:astro- 
ph/0509531]; J. Y. Shen, B. Wang, E. Abdalla and R. K. Su, Phys. Lett. B 609, 200 (2005) [arXiv:hep- 
th/0412227]; Z. L. Yi and T. J. Zhang, Mod. Phys. Lett. A 22, 41 (2007) [arXiv:astro-ph/0605596]; Y. Z. Ma, 
Y. Gong and X. Chen, Eur. Phys. J. C 60, 303 (2009) [arXiv:0711.1641 [astro-ph]]; Q. Wu, Y. Gong, A. Wang 
and J. S. Alcaniz, Phys. Lett. B 659, 34 (2008) [arXiv:0705.1006 [astro-ph]]; M. Li, X. D. Li, S. Wang, 
Y. Wang and X. Zhang, JCAP 0912, 014 (2009) [arXiv:0910.3855 [astro-ph.CO]]; M. Li, X. D. Li and X. Zhang, 
arXiv:0912.3988 [astro-ph.CO]. 

[19] H. Wei and R. G. Cai, Phys. Lett. B 663, 1 (2008) [arXiv:0708.1894]; I. P Neupane, Phys. Rev. D 76, 
123006 (2007) [arXiv:0709.3096 [hep-th]]; J. Zhang, X. Zhang and H. Liu, Eur Phys. J. C 54, 303 (2008) 
[arXiv:080 1.2809 [astro-ph]]; Y. W. Kim, H. W. Lee, Y. S. Myung and M. I. Park, Mod. Phys. Lett. A 



14 



23, 3049 (2008) [arXiv:0803.0574 [gr-qc]]; J. P. Wu, D. Z. Ma and Y. Ling, Phys. Lett. B 663, 152 

(2008) [arXiv:0805.0546 [hep-th]]; J. Cui, L. Zhang, J. Zhang and X. Zhang, Chin. Phys. B 19, 019802 
(2010) [arXiv;0902.0716 [astro-ph.CO]]. L. Zhang, J. Cui, J. Zhang and X. Zhang, Int. J. Mod. Phys. D 
19, 21 (2010) [arXiv:09 11.2838 [astro-ph.CO]]; X. L. Liu and X. Zhang, Commun. Theor. Phys. 52, 761 

(2009) [arXiv:0909.4911 [astro-ph.CO]]; X. L. Liu, J. Zhang and X. Zhang, Phys. Lett. B 689, 139 (2010) 
[arXiv: 1005.2466 [gr-qc]]. 

[20] X. Zhang, Phys. Rev. D 79, 103509 (2009) [arXiv:0901.2262 [astro-ph.CO]]. 

[21] C. J. Feng and X. Zhang, Phys. Lett. B 680, 399 (2009) [arXiv:0904.0045 [gr-qc]]; C. J. Feng, arXiv:0806.0673 
[hep-th]; C. J. Feng, Phys. Lett. B 670, 231 (2008) [arXiv:0809.2502 [hep-th]]; L. Xu, W. Li and J. Lu, 
Mod. Phys. Lett. A 24, 1355 (2009) [arXiv:0810.4730 [astro-ph]]; C. J. Feng, Phys. Lett. B 676, 168 (2009) 
[arXiv:08 12.2067 [hep-th]]; K. Y. Kim, H. W. Lee and Y. S. Myung, arXiv:08 12.4098 [gr-qc]. 

[22] R. G. Cai, Phys. Lett. B 657, 228 (2007) [arXiv:0707.4049 [hep-th]]. 

[23] R. G. Cai, B. Hu and Y. Zhang, Commun. Theor. Phys. 51, 954 (2009) [arXiv:0812.4504 [hep-th]]. 
[24] L. Pasquini et al, Proc. Int. Astron. Union 1, 193 (2005). 

[25] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 180, 330 (2009) [arXiv:0803.0547 [astro-ph]]. 



